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Abstract 

375.0ptWe investigate the local energy decay of solutions of the Dirac equation 
in the non-extreme Kerr-Newman metric. First, we write the Dirac equation as 
a Cauchy problem and define the Dirac operator. It is shown that the Dirac 
operator is selfadjoint in a suitable Hilbert space. With the RAGE theorem, we 
show that for each particle its energy located in any compact region outside of 
the event horizon of the Kerr-Newman black hole decays in the time mean. 

1 Introduction 

The Kerr-Newman metric is the most general stationary solution of Einstein's field 
equations and has the physical interpretation of a massive charged rotating black hole. 
The physical parameters M, Q and a = J/M have the interpretation as mass, electric 
charge and angular momentum parameter. 

As in flat spacetime, a spin-i particle is described by a Dirac equation. In the Kerr- 
Newman metric, the Dirac equation is given by the coupled system of partial differential 
equations (see, e.g., [Pag76], [Cha98]) 

($H + 2l)§ = (1) 

where 



m := 



/ 


imr 







° \ 







—imr 














—imr 





\ 










imr ) 



21 



/-» 








£*/\ 





s 










£*/ 


-S 













3) / 



* winklmeierCSmath. unibe. ch 
t yamadaos@se. ritsumei .ac.jp 



(2) 



1 



and 



D := 



am cos 6, 
dr A 



d 



d 



d_ 

89 



cote 



(r 2 + a 2 ) i— + a i— + eQr 
. <9 



a* 



9^ 
1 



(3) 

on(r+,oo), (4) 



asmf l 



d 



dt sin dip 



on (0, 7r). 



(5) 



The physical parameters m and e are the mass and the electric charge of the spin-^- 
particle which is described by the wave function $ . 
The functions A and £ are given by 

A(r) = r 2 - 2Mr + a 2 + Q 2 = (r - M) 2 + a 2 + Q 2 - M 2 , 



a 2 cos 2 0. 



According to the black hole parameters M, Q and a, three cases can arise: If a 2 + 
Q 2 — M 2 < 0, then the function A has exactly two distinct zeros; this case is called 
the non-extreme Kerr-Newman metric. If a 2 + Q 2 — M 2 = 0, then the function A 
has exactly one zero; this case is referred to as the extreme Kerr-Newman metric. In 
both cases the metric is interpreted physically as the spacetime generated by a charged 
rotating massive black hole. If a 2 + Q 2 — M 2 > 0, then the function A has no zeros; in 
this case the Kerr-Newman metric has no interpretation as the metric of a black hole. 

It is a remarkable fact that the Dirac equation (1) can be separated into ordinary 
differential equations. This was first shown by Chandrasekhar ([Cha76]).The first step 
is to separate the dependence on the coordinates t and <p by the ansatz 



(6) 



uj G M. We call w an energy eigenvalue of equation (1) if it has a solution of the 
form (6) which is square integrable in the sense explained below in section 2. 
The separation of the Dirac equation into ordinary differential equations is achieved 
by the ansatz 



*(r, e,<p,t) 



/X_(r)S + (6)\ 
X+(r)S_(0) 
X + (r)S+(e) 

\X-(r)S-(9)J 



with k £ r L+\. One obtains two differential equations, the so-called radial equation for 
the radial coordinate r and the so-called angular equation for the angular coordinate 
6. Both equations can be written as eigenvalue problems for selfadjoint operators 
in appropriate Hilbert spaces. The spectrum of the corresponding angular operator 
consists of simple and discrete eigenvalues which are unbounded from below and above 
(see, e.g., [Win06]). Bclgiorno and Martellini [BM99] showed that the spectrum of 
the radial equation comprises all of the real axis. Recently, it was shown that in 
the non-extreme Kerr-Newman metric no embedded eigenvalues exist. In the extreme 
Kerr-Newman case, however, embedded eigenvalues can exist (see [Sch04] and [WY06]). 

In this paper we investigate the local energy decay of solutions of the Dirac equation (1). 
We assume that the non-extreme case holds, that is, there arc no energy eigenvalues. 
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We do not employ the separation ansatz described above, but we rewrite the Dirac 
equation as a Cauchy problem. After some transformations, we obtain the equation 



i—V{x,6,tp,t) 



(7) 



where S(x, 9) is a bounded and boundedly invertible 4 x 4-matrix and Sj is the formal 
Dirac operator associated with the Dirac equation (1). We show that the expression Sj 
on the right hand side has a selfadjoint realisation H in an appropriate Hilbert space 
(Theorem 2.1). The operator H can be written as the sum over partial wave operators 
H( K \ k G Z. In order to apply the RAGE theorem to the Dirac operator, we show in 
Theorem 3.3 that it enjoys a Rellich type property. One of our main results in this 
paper is Theorem 4.2 which shows the local energy decay of the wave functions in 
the time mean: Let be the group associated with the skew-selfadjoint operator 
—iS" 1 !!^'. Then, any function / which is subjected to the time evolution given by 
(7) satisfies 



lim — 

T-oo 2T 



R, r-K 



R JO 



U^(t)f(x,6) 



dx d6 



dt = 0. 



for every R > 0. Since the expression in the square brackets is related to the energy of 
the particle described by / at time t in the region Q := (—R,R) x (0, if), this result 
is related to the decay of the particle's energy in f2o- The above result remains valid 
if the bounded electric potential r t+a 2 given by the background metric is substituted 
by a possibly unbounded potential q satisfying 

q G C 1 ((r + , oo)), lim ^ exists, q'(r)=O0r), as r — > oo. 

r — >oo t 



The local energy decay has also been investigated by Finster, Kamran, Smoller and 
Yau ([FKSY03]). They restrict the Dirac operator to an annulus outside the event 
horizon and impose Dirichlet type boundary conditions so that the restricted operator 
has purely discrete spectrum. Hence the propagator associated with it is a sum of 
projections. Extending the radii of the annulus to the event horizon and infinity, 
respectively, the sum of the propagator representation becomes an integral. 
Our approach, however, deals directly with the Dirac operator in the exterior of the 
black hole. Using the fact that the partial wave operators H^ K > have no point spec- 
trum, the local energy decay of wave functions with initially compact support follows 
by methods from scattering theory. In particular, we do not impose any additional 
boundary conditions. 

2 Transformation of the Dirac equation to a 
Cauchy problem 

To investigate the time evolution of solutions of the time-dependent Dirac equation we 
rewrite it as a Cauchy problem. 

In the following, we use the Pauli matrices to abbreviate notation: 




3 



In addition, let 

*-(; !)■ "-(« -".)■ :)--^ 

For the (^-coordinate we always assume the boundary condition ^(r, 0, <p) = \&(r, 0, 95 + 
2tt). 

In the rest of the paper we use the notation 

Q 3 := (-00,00) x (0,7r) x [0,2tt), 
2 := (—00,00) x (0,7r). 

Theorem 2.1. Assume that the function A has at least one real zero (that is, either 
the extreme or the non- extreme Kerr-Newman case holds) and denote the largest zero 
by r+. Define the new radial coordinate x G (—00, 00) by 

dx r 2 + a 2 



dr A(r) ' ^ 

To keep notation simple, we often write r instead r(x). 

A function ^ is a solution of the Dirac equation (1) if and only if 

i^* = 5- 1 U* (9) 

where 

*(x,9,<p,t) = (smfl)- 1 / 2 *^),^), (10) 



S 



aVAsin^ / a-? 

1+—, 7T~ 



-0-2 



55 = Sji +$j 



2- 



\ — <9r r 2 + a 2 \ 90 "sinfl 

^ = -^T^C 2 / 2 )( ai ^ + cgr 

mr\f~K ( 7 2 \ am\f~Kcos6 f i/ 2 
r 2 + a 2 / r 2 + a 2 y— i/2 

Moreover, the operator Hq defined by 

= D(ff ) = ^T^a) 4 (11) 

is essentially selfadjoint in the Hilbert space Jzf 2 (il3) 4 and the operator 

S _1 flo 

is essentially selfadjoint in the Hilbert space 

J^j(ft 3 ) 4 :=J^ 2 (ft 3 ) 4 wit/i scalar product ( • , • ) s 

given by 

r /*27r 

(*, *) s = (*, 5$) = / / / (*(x,0,tp), S(x,9) *{x,6,<p))& dx dOdip; 
Jrjo Jo 



accordingly, the norm on ^^(^3) will be denoted by 



s- 
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Let H be the closure of Hq. Then S 1 H is the closure of S 1 Hq 1 and we call S 1 H 
the time-independent Dirac operator in the Kerr-Newman metric. 

Remark 2.2. Observe that for each (x, 6) € (—00,00) x (0,7r) the matrix S(x,9) is 
bounded and boundedly invertiblc since 



a\fK sin 9 



r(xy 



< 



\a\r(x) 



r(x) 2 



< 



1 



Therefore, the norms on .i? 2 (r2 3 ) 4 and Jz?|(Sl 3 ) 4 arc equivalent. 

Proof of Theorem 2.1. First we rearrange equation (1) such that all time derivatives 
are on the left hand side and all other terms are on the right hand side, thus we obtain 



-10-3 
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Multiplication from the left by ( ; 



avAsinf? /0? 



-0-2 



-03 







yields 



i—^(x,e,tp,t) 

at 



A 



0-3J dx r 2 + a 2 



d_ 

<Ji J *~d9 



(r 2 + a 2 ) sin( 



-02 



d 



o~2j dip r 2 + a 2 



h) \ a % +CQr 



mr\f~K 



1 2 U 



amy A cos 9 ( i J 2 

-i/ 2 



- (j5i + J5 2 )*(a:,0,<M), 
thus the first assertion of the theorem is proved. 

By Remark 2.2 the norms on Jz? 2 (0 3 ) 4 and Jz?J(f2 3 ) 4 are equivalent. Hence, TJo is 
essentially sclfadjoint in _Sf 2 (H3) 4 if and only if iS _1 i?o is so in ^fj(il 3 ) 4 . Therefore, it 
suffices to show that Ho is essentially sclfadjoint. 

Since ^ with the boundary condition stated at the beginning of this section has the 
complete system of eigenfunctions {e~ 1KLp ) K e"L+i/2 and all the differential expressions 
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in the equation above commute with each ^ e «5f 2 (SI3 x M) has a representation 



= E e- iK ^(x,e,t). 

K6Z+1/2 

\I/ satisfies (9) if and only if each 1 i> < - K > satisfies 

i-^-^W = (12) 
dt K ' 

where fi( K \ Sy^ and ^2 are obtained from Sj, Sji and Sj 2 by substituting 

^~d<p with 

Moreover, Hq is essentially sclfadjoint if and only if each partial wave operator 
(the restriction of H to the partial wave space represented by e~ 1Kip ) is essentially 
sclfadjoint. 

In Proposition 2.4 below, it is shown that h[ K Q, defined by 

= f)i K) * (K) , = ^{VL 2 )\ (13) 

is essentially selfadjoint in Jz? 2 (r22) 2 with the usual scalar product. 
Since both f)i and 6> have diagonal structure and Sj 2 defines a bounded sclfadjoint 
operator on _SfJ(f22) 4 , the assertion on the essential selfadjointness of Hq , k£ Z, and 
therefore of iJo follows. □ 

Lemma 2.3. For k 6 Z + 1/2 fe£ 

A a k /o 1 



and define the unitary matrix 

W := 



-1 0/ 96» sinfl U 



f 

1 



The operator A K is the angular part of the Dirac equation (1) arising from Chan- 
drasekhar's separation process in the case a = 0. 

(i) For each n £ Z + 1/2, the operator A K with domain "^"((O, it)) 2 is essentially 
selfadjoint in J;? 2 ((0, 7r), dfl) 2 . FKe denote its closure again by A K . It is compactly 
invertible and its spectrum consists of simple eigenvalues only, given by 

Kn ~ sign(m) - \ + \m\) , m = ±1, ±2, • • • . 

We denote the corresponding normalized eigenf unctions by g 1 ^ . 

(ii) For every m £ Z\ {0} we have \1 m = — A£, and 03.0^ = —g1 m . 
(hi) WVUW" 1 = W-M K W = i(r 3 ^- re - 

9\ y/A /.d\ VA 
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Proof. For (i), we refer to [Win06]. Statements (iii) and (iv) are easily verified by 
direct computation. To prove (ii) we observe that all eigenvalues of A K are simple and 
that 



1 Ul 
-lMo -1 



1 

-i 



1 

-i 



9m - -\n a 39 m - — A -m (J 35 m - 



□ 



Proposition 2.4. The minimal operator h[ K Q defined in (13) by 



. d 
-0-31 TT 



n2 



=^°° ((-00,0c) x (0,tt)) 2 
is essentially selfadjoint in the space J£ 2 (Vl2) ■ 

Proof. Obviously f)^ is hermitian and is symmetric. Hence, the lemma is proved 
if we have shown that ker(/ij q ± i) = {0}. The equality h^^ip = ±iip is equivalent 



to 



with the matrix VT as defined in Lemma 2.3. 

If we set (p := W~ lr ip, then, by Lemma 2.3, the equation above is equivalent to 



03 -1 



d_ 

dx 



H + a 1 



(14) 



Since the set of eigenfunctions of the operator A K , denoted by g% as in Lemma 2.3, is 
complete in jSf 2 ((0, 7r), d8) 2 , we can expand the function ip in (14) as 

m£Z\{0} 

with functions £ ro e ^f 2 ((— 00, 00), dx). Thus (14) yields 
±i E £. m (x)g m (0) 



mGZ\{0} 



mez\{o} 
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Taking the scalar product with g^iO) in C 2 and integrating with respect to 8 yields on 
the left hand side 

±i r E U(9^(o)9^e)) c2 de = ±i E U f (g h m {o)9^)) c ,de 

J ° meZ\{0) mSZ\{0} J 

= ±i£ M 

and on the right hand side 

E iimWff- m (^) + M(m(^(fl), 9^)) dB 
OX T ~\~ CL j ^2 

E '^rnix) J (g- m (6), g^0))c d6 



m6Z\{0} 



E Ms) f\g m {6), gM)& dd f\g m {9l git (9))c> d9 



Ox r z + a 1 



Hence we obtain 



±i^W = i^W + ^fpW- (15) 



An analogous procedure with <7_^ gives 



±i^(x) = + (16) 

ox r z + or 



Combining equations (15) and (16) we obtain 





(17) 



The multiplication operator is bounded on Jzf 2 ((— oo, oo), dx) and it is well known 

that the operator ( . ^_ 1 J** "] with domain <r ^ x) (— o3,co) 2 is essentially selfadjoint in 
Jzf 2 ((— oo, oo), dx) 2 , see, for instance, [Wei87, Ch.6.B)] . Hence the same is true for 



A2 n/A 



/ r-t+a-z A f^ /S . J and equation (17) has only the trivial solution in Jz? 2 ((— oo, oo), dx). 
' 35 ~ ; 'i +a 'i ' 

This, however, implies that also equation (14) has only the trivial solution and the 
lemma is proved. □ 

Let be the unique selfadjoint extension of Sj^ (see (12)). The following result is 
proved in [Sch04]. 

Proposition 2.5. For each half-integer k, the spectrum a(H^) coincides with R, and 
has no eigenvalues. 
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3 Rellich property 



Definition 3.1. Let f2 be an open set in 1™ and T a selfadjoint operator on Jz? 2 (£l). 
We say that T has the Rellich property, if the following is true: Let F C T>(T) be such 
that there exists a positive constant K such that 

ll/ll 2 + !|77|| 2 < K, feF. 

Then F is precompact in L 2 oc (Q!) for every bounded open set fi' C fi, that is, every 
sequence (/„)„ C F has a subsequence [fn k )k which converges strongly in j£? 2 (f2'). (As 
usual, we identify functions f € F with the corresponding functions /|q< € Jzf 2 (f2').) 

The classical result that differential operators in Jz? 2 (R d ) have the Rellich property is 
cited in the appendix. 

The aim of this section is to show that Sj^ has the Rellich property. First we proof 
the following technical lemma needed for Theorem 3.3. 

Lemma 3.2. Let W and A K , n € Z+ 1/2 be defined as in Lemma 2.3 and let be 
the Dirac operator in the n-th partial wave space (see Theorem 2.1). 
Then there exists a positive constant K such that 



2 

E 

3=1 



d6 dx 



(r 2 + a 2 ) 2 

< K(\\S- l H^f s 

for all * = *(*i, * 2 ) G ViS^H^)) where *j £ «^j(^ 2 ) 2 7 j = 1, 2. 
Proof. First assume e ^ oo (fi 2 ) 4 . Observe that 



1*111) (18) 



< 



a sin^ \/A 



< 2, 



(r,0) e (-00,00) x (0, 7r). 



With the notation of Section 2 we write 5 1 i^( K ) as the sum of an unbounded operator 
anc j a b oun d ec j 5 _1 7J 2 K ' ) and obtain 

Iis-^^H 2 > 1 ii.s- 1 ^*!! 2 -n.s- 1 ^^!! 2 
1 



i] 



S~ 1/2 H[ K) ^{x,9)\\ 2 cl dxde - ws^H^hwl 



Recall that 



(«) _ / ^3 



1 

> - 

- 4 



c 3 



H l f ) ^!(x,e)\\l i dxde - HS" 1 ^!! 2 
\/A (w~ x A K w 



dx J r 2 + a 2 



-W~ A K W 



(19) 



0-3 



-03 



9a; 



VA /i^ 



-i A, 



by Lemma 2.3. Since the first matrix in the above representation is unitary, we obtain 
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\H[ K) y(x,9)\\l 2 dxd6 



= E 

3=1 
2 

= E 



• 9 T/ 



• 9 T , 



a" 



*(x,6>) 



dd dx 



'A 

- T —^iA K * J {x,< 



d6 dx 



9 T 



A 



c 2 



(r 2 + a 2 ) 5 



| A i^, 6011c 2 



ciflcfa: (20) 



c 2 



where in the last step we have used that — i d x is symmetric and A K commutes with 
— \d x and performed integration by parts. There are no boundary terms since \I/ has 
compact support. 

There exists a positive constant C such that 



A 



dx r 2 + a 2 



dr d 



dx dr r 2 + a 2 
A 



A 



r 1 + a 2 



dr r 2 + a 2 



r 2 + a 2 



2rA 



r - M 



(r 2 + a 2 ) 2 r 2 + a 2 



< C- 



(21) 



since the term in | ■ | is continuous and tends to zero for x — > ±oo. Moreover, for every 
£ > 



d 



Jl \ dxr 2 + a 2 



A K *i 



> -II* 



Jiiif 2 (n 2 ) 2 



dx r 2 + a 2 j K 



_e 2 2 1 
- 2 " J '"^ 2 (^2) 2 2 £ 2 



.S? 2 (0 2 ) 2 



(22) 



3f 2 (a 2 ) 2 



Hence, from (19), (20), (21) and (22) we obtain 



J 11*11^(0, 



(r 2 + a 2 ) 2 



(1 ^^"^fe 9)112, 



eft? e?a; 
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If we choose e small enough such that < 1 and observe that 1 1 ^ 1 1 (o 2 ) 4 — 
H^-!/ 2 ]] 2 ||\I>|| S then we can choose an K > large enough such that the estimate 
in the assertion holds. 

For any "J" £ V(H^) there exists a sequence ( X P„)„ such that 

in _S?5 (02)- Notice that for this sequence also the left hand side of (18) converges by the 
dominated convergence theorem. Therefore, assertion holds for all \I> £ V{H^). □ 

Note that K does not depend on \& but only on the J^-bound C of \fK-^ 



dr r + a 2 

Theorem 3.3 (Rellich property). Let k 6e any feaZ/ integer, R > and sei 

n 2 ,j l :=(-i2,ii)x(-7r > 7r), (23) 
Let (*„)„ C V(H^) such that 

for some constant Kq > 0. Then there exists a subsequence (^ ni )i and a $ £ 
^ 2 {^2,r) A such that 

&n t — > $ as £ — > 00 

zn^j^) 4 . 

Remark 3.4. As usual, we identify elements \& £ „Sf 2 (f2 2 ) 4 with elements of 3/ £ 
-S? 2 (^2._r) 4 by restriction. The assertion of the Theorem implies that there exists a 
subsequence (^ ne )e and a $ £ Jz? 2 (fi 2 ) 4 such that XR®n t — > Xfl.^ as ^ ~ * 00 m 

if|(n 2 ) 4 . 

Proof of Theorem 3.3. Since the norms on Jz?J(f2 2 ) 4 and Jz? 2 (0 2 ) 4 are equivalent, it 
suffices to consider convergence in the latter space. Let >J>„ =: (tp^, tp 2 ) with ip^ £ 
«5f 2 (il 2 ) 2 , j = 1, 2. Lemma 3.2 and the existence of a positive constant 5 > such 
that 

, 2 ^ 2 n 2 > 8, xe (-R, R) 
(r 2 + a 2 ) 2 

imply that there is an K > such that 



2 /.« 

E, 



d 



+ \Mi(x,9)\\' c2 dxde 

c 2 I 



< K(\\* n \\ 2 + \\H^* n \\ 2 ) < KKl (24) 

for n £ N. Since Jzf 2 ((0,7r), dO) 2 has a basis of orthonormal eigenfunctions (jJ5,) m of 
with corresponding eigenvalues AJ^, m £ Z\ {0}, see Lemma 2.3. Hence ipn can be 
expanded in the double Fourier series 



2i? J-RjQ 

/-7T 



'oils- 



v. in 
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Moreover, inequality (24) yields 

£ [(vn/R) 2 + (\ h m ) 2 } \a v , m , n \ 2 < KK . 

For fixed (is, m) the sequence (a v ,m,n)n is a bounded sequence, hence it contains a 
subsequence such that (a v , m ,ni)m is a Cauchy sequence for any (v, to) by a diagonal 
series argument. Thus for arbitray L > it follows that 



2R 



R fK 



R JO 



l^-^/dxde 



^ l-R. /-7T 

7775 / / 1^4 - V'i. | 2 dxdfl = V la^m.n, - a^, m ,„ 3 

^RJ-RJo fZ 



E 



i/.m 

2 



+ E 

(^7T/i?) 2 +(A^) 2 >L 

< E i 

(i/7i-/fl) 2 + (A^) 2 <L 



[(fTf/i?) 2 + (A^) 2 ] \a v , m ,n e - Uv,m, nj f 

(w/Ry + (a« y 

1 + L' 



which shows that (ip ne )e is a Cauchy sequence in Jz? 2 (0 2 ) 2 . 

With the same argument we find that also the sequence (tp 2 t )£, contains a convergent 
subsequence in J2f 2 (f2 2 ) 2 - □ 



4 Weak local energy decay 

Since H^ K \ k £ Z \ {0} is selfadjoint in «5fJ(02) 4 , it is generator of a unitary group 

t/ (K >(t) :=exp(-it5- 1 J ff (K) ). (25) 

Theorem 4.1. Let K G Z\{0}, i? > and * G V(H^). For any sequence (t n )„ C R 
£/iere e:risi a subsequence (t ni )l and <E> G J2? 2 (f22) 4 such that 

t/ (K) (i«J*o^$ as 



oo 



m ^ 2 (0 2 ) 4 . 

Proof. Let := U^(t n )^ . Then, for all n G N, 

||¥„|| S + ll^- 1 ^)* n || S = ||tfW(t„)*o||s 

= ||^W(* n )*o||s 

= ||* ||5 + ||«S- 1 i/ (K) *oi|5. 

Thus the assertion follows by Theorem 3.3 with Kq = \\^o\\ + H-ff^Vto 



C/W(i n )5- 1 ffW* ||s 



□ 



Theorem 4.2. For every R > and $ G .^(f^) time mean of the localization of 
$ m (— i2, i?) x (0, 7r) is to zero: 



lim — 

T^oo 2T J_ T 



R fK 



-R JO 



u {K) (t)$>( x ,e) 



dx dO 



c 4 



dt = 0. 



12 



Proof. It is well-known that in the non-extreme Kerr-Newman case the operator 
has no eigenvalues ([Sch04], [WY06], see also Proposition 2.5). Let xr( x ) be the 
characteristic function of (—R,R). Then, the Rellich property Theorem 3.3 implies 
that xn(H^ + i)^ 1 is a compact operator. Therefore, the assertion follows from the 
RAGE theorem, see [RS79]. 



□ 



The decay of the partial waves implies the decay of the solutions of the original Cauchy 
problem (9): 

Theorem 4.3. Let * € Jz?j(ft 3 ) 4 . Its Fourier expansion 

«ez+(i/2) 

converges strongly in jSf^ (O3) 4 and 



Then 



||*«(M)Hl < 00. 

«ez+(i/2) 



reeZ+(l/2) 



zs £/ie unique weak solution of 

. 8 



dt 



f(0) = * 



(c/. (9)) m Jz?J(f23) 4 and satisfies 



lim — — 



it /.it /-27T 



RJO JO 



dx d6 dip 



dt = 



/or any R > 0. 

Remark 4.4. The above results remain still valid if the electric potential eQr in (4) is 

substituted by any real- valued C -function q(r) defined on [r + , 00) such that lim 

) — >oo r 

exists and q'(r) = O(l) as r — > 00. 

A Equivalent representations of the 
time-independent Dirac operator 

Remark A.l. If we do not apply the transformations (8) and (10) to the wave function 
the time- independent Dirac equation takes the form 



at 



(26) 



with f) = f}\ +£)2 with f)2 as in Theorem 2.1 and 
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r A . d VA / ./ d cot0\ 1.5 
f)i = -03^— — ni t; 5— — o o-i 1 ™ H ^— + 0-2^^71 



2 / sin dip J 
$j is formally symmetric in the Hilbert space 

i|:=J^ 2 ((r + ,co)x(0,7r)x[0,27r)) 

with the scalar product 

POO P 7T p27T A 

/ / ($(r,e, V ),5(r,fl,vj)$(r,fl,vj)>c«-2X- ^dr sind d9 dtp. 

Jr + Jo JO r " + a 

Remark A. 2. In Theorem 2.1 we have equipped the Hilbert space Jz? 2 ((— 00,00) x 
(0, 7r) x (0, 2ir)) with a scalar product ( • , such that is formally symmetric. Al- 
ternatively, we can work with the usual scalar product, but then we have to define the 
Dirac operator in the following way: 

H V = S- 1/2 SjS~ 1/2 ^, V(H ) = <r o °°((-oo, 00) x (0, tt) x [0, 2tt)) 4 

on Jz? 2 ((— 00, 00) x (0, 7r) x [0, 7r)) 4 equipped with the usual scalar product 

(*, $r= / / / (*M.¥>), 3>{x,e,<p)) c *dxd0dip. 

JR JO Jo 

Note that <S leaves the space ^ U X3 ((— 00, 00) x (0, 7r) x [0, 27r)) 4 invariant. 

B Rellich's theorem 

Theorem B.l (Rellich's theorem). Let d G N. A' C M d compact and a bounded 
open neighbourhood of K . Further let s > t G N. Then for each c > the set 

F :={ue H s (R d ) : \\u\\ H s (Rd) < c, supp(u) C A} 

is precompact in iJ*(0), i.e., /or eac/i sequence (f n )n C A t/iere exists an /o £ if*(R ) 
and a subsequence (f nk )k such that f nk — > /o m H (Q,). 

Proof. See, e.g., [Miz73]. □ 
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